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SCATTERING  BY  WATER  WAVES  GENERATED 
BY  A  MOVING  PRESSURE  POINT 

R.  H.  Ott,  G.  A.  Hufford  and  A.  Q.  Howard,  Jr. 


ABSTRACT 


,  An  analysis  for  scattering  of  scalar  waves  in 
air  by  a  wake  on  the  surface  of  water  is  developed.  The 
wake  on  the  water  surface  is  generated  as  a  concentrated 
pressure  point  travels  over  the  free  water  surface.  The 
wake  has  characteristics  resembling  a  diffraction  grating. 

The  grating  produces  maximum  constructive  interference 
(resonance)  in  the  backscattered  wave  under  certain 
conditions.  At  resonance  the  acoustic  wavelength,  water 
wavelength  and  angle  of  incident  satisfy  a  simple  geo¬ 
metrical  condition. 

AVi 

1.  INTRODUCTION 

An  analysis  for  studying  the  scatter  of  scalar  waves  in  air  by 
a  wake  on  the  surface  of  water  is  developed.  The  wake  on  the  water 
surface  is  generated  by  a  cone  nitrated  pressure  point  traveling  with 
uniform  velocity  over  deep  water.  The  only  surface  force  i,  pressure, 
the  only  body  force  is  gravity  and  the  motion  is  irr otational. 

A  number  of  authors  have  investigated  the  wake  left  by  a  pres¬ 
sure  point  as  it  travels  over  the  water.  The  wake  is  largely  confined 
to  the  sector  bounded  by  two  cuspidal  lines  which  form  angles  of  sin"1  (£) 
with  the  track  of  the  pressure  point.  The  disturbance  drops  off  expon¬ 
entially  outside  the  cuspidal  lines.  This  angle  was  first  determined 


by  Lord  Kelvin  (1887).  Within  this  sector  the  wake  consists  of  a  system 
of  waves  that  envelope  the  hull  lengthwise  and  is  interwoven  with  a 
system  of  transverse  waves.  The  amplitude  of  the  surface  within  this 
sector  has  been  treated  by  a  number  of  authors  (e.g.  ,  Peters  (1949). 
Ursell  (I960),  Hufford  (1972)). 

In  this  paper  we  introduce  the  boundary  value  problem  for  the 

velocity  potential  5(x,y.z)  and  the  surface  elevation  C  <*, y) .  The 
velocity  potential  satisfies  Laplace's  equation.  The  velocity  potential 
and  surface  elevation  together  satisfy  two  boundary  conditions:  1)  the 

pressure  equation  and  2)  the  equation  of  continuity  across  the  free 
surface. 

Sommerfeld  (1950)  is  a  good  introduction  to  the  theory  of  waves 
and  the  equations  of  hydrodynamics  applied  to  our  problem. 

The  approach  to  the  solution  to  Laplace' s  equation  plus  boundary 
conditions  is  through  the  two-dimensional  Fourier  transform  (Havelock 
(1918),  Hogner  (1923),  Peters  (1949)).  These  equations  are  solved  in 
terms  of  the  transformed  potential  *(S,n.*)  and  surface  elevation 
Z(§,ri).  The  final  step  consists  of  identifying  the  transform  Z  (§,q) 

as  the  solution  to  our  scattering  problem.  A  simple  expression  for  the 
scattered  field  is  obtained  which  depends  upon  the  acoustic  wave 
number  k,  the  water  wave  number  kw  and  the  direction  of  the  source 
and  observer  relative  to  the  track  of  the  pressure  point.  This  expres¬ 
sion  has  characteristics  resembling  the  field  scattered  by  a  diffraction 
grating.  The  interference  phenomena  in  this  problem  are  explained 
physically  by  reflections  from  neighbouring  curves  having  water  wave 
phase  differences  which  are  multiples  of  2tt. 
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2.  ANALYSIS 


Consider  a  monochromatic  scalar  wave  having  time 

dependence  incident  on  the  surface  z  =  £(x,y)  in  figure  1.  Assume  the 
potential  C p  satisfies  a  sound-hard  boundary  condition  (corresponding  to 
a  vertically  polarized  incident  wave  in  the  electromagnetic  problem) 


acp 

dn 


=  0  , 


on  i 


(1) 


where  differentiation  is  with  respect  to  the  outward  normal.  The 
scattered  potential  from  Helmholtz's  formula  is 


-/<P)  ■  L  J  ,(Q)  |2  da 

C 


(2) 


where  P  is  the  location  of  the  source  and  observer  for  the  case  of 
backscatter  and  Q  is  the  integration  point  and  £  is  the  surface  of 
integration.  The  free-space  Green's  function  is 


-ikr~ 
e  ® 

c  -  —  <3> 

with 

r2  =  PQ  .  (4) 

Neglecting  amplitude  terms  of  the  order  (kr2)-1  ,  we  obtain  from  (3) 


dG 

dn 


-ikG 


(5) 


We  assume  that  near  grazing  the  field  at  the  integration  point 
Q  is  the  unperturbed  field  of  a  point  source  at  P  radiating  over  a 
perfectly  conducting  flat  surface;  i.e.  , 
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(6) 


•  -  ii\t2 

cp(Q)  =  2cp  (Q)  =  2  - - 

r3 

The  elemental  area  da  projected  onto  the  x-y  plane  becomes 


e  •  e 
— n  —  z 


where  e^  is  the  outward  directed  normal  to  Q  in  figure  1  and  e  is 
a  unit  vector  along  the  positive  z -axis.  The  normal  derivative  of  r2 
in  (5)  is 

d  ra 

an  =  —  r2  *  — n  (8) 

Wlth  £rs  a  unit  vector  from  the  observation  point  P  to  the  integration 
point  Q.  The  unit  normal  is 


-n=(-TC  +  Sz)l 


V  C  -  ,  d£ 

-Tc  "  -X  bZ  +  £y  a7  (10) 

and  V^C  is  shown  in  plan  view  in  Fig.  1.  That  is,  is  the  projec¬ 

tion  (to  within  the  normalization  •  e  J  of  £n  onto  the  x-y  plane.  If 
the  surface  elevation  £  is  small  compared  with  the  acoustic  wavelength 
X  and  the  slopes  dC/dx  and  d£/dy  are  small  compared  to  unity  we 
find  for  an  observer  close  to  the  x-y  plane 

(11) 


with  ek  2:  a  unit  vector  from  P  to  the  origin  in  figure  1.  From  (8), 

(9),  (10)  and  (11) 


I?-  “  (*k  ‘  V)  I  4» 
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Substituting  (3),  (5),  (6),  (7)  and  (12)  into  (2)  gh 


-i2kr- 


— 2  dx  dy 

r2 


*  (p)  =  TtT  J^k‘  -Tcj  7F  ixdy 

where  C  is  the  projection  of  Q  onto  the  x-y  plane. 


F rom  figure  1 . 


r2  -  R  +  e ,  •  r 
— k  — 


r2  "  R  +  (x  cosfl  -  y  sin 9) 


Substituting  (15)  into  (13)  and  making  the  usual  approximation  regarding 
the  amplitude  factor  r2_1  we  obtain 


9  (P)  = 


. ,  -i2kR  .  v 

J  (*k  ■  -tc) 


-i2k(x  cos0  -v  s  in© ) 
e  y  'dxdy  (16) 


which  is  a  two-dimensional  Fourier  transform  of  f  efc  •  vTCj  .  The 
surface  is  assumed  to  extend  from  -®  to  + 00  along  the' x-  and 


y-axes. 


F  rom 


— k  *  — =  COS0  i  - 


and  the  property  for  the  Fourier  transform  of  the  derivative  of  a 
function,  we  obtain 


,„s.  -  kz  -i2kR 

*  (p)  =  PpJ  e  Z(?  .  n  ) 


where 
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Z(S,T])  =  J  dx  j  dyC(x.y)  e‘l(§x  +  T1y)  (19) 

-00  -  00 

is  the  Fourier  transform  of  the  surface  elevation  G(x,y).  We  have 
assumed  in  (16)  that  Q*  is  infinite  in  extent  and  at  the  same  time  to 
arrive  at  (16)  we  make  approximations  regarding  the  phase  kr2  in 
(14)  and  the  amplitude  factor  rg"1  .  These  approximations  can  be  made 
rigorous  using  the  theory  of  generalized  functions  (Hufford,  1972). 
Hufford  shows,  by  considering  a  sequence  of  surfaces  {Cn}  whose 
limit  is  C,  that  in  the  limit  as  G^  ~*co  the  scattered  field  approaches 
the  result  in  (16). 

The  variables  §  and  T)  in  (19)  become 

5=ZkCOs9  (20) 

T]  =  -2  k  sin  6 

when  they  are  used  to  denote  the  components  of  the  acoustic  wave 
number;  otherwise,  they  will  be  thought  of  as  the  components  of  the 
water  wave  number. 

The  result  in  (18)  shows  that  finding  the  field  scattered  by  a 
rather  complicated  surface  G  is  reduced  to  the  problem  of  finding  the 
Fourier  transform  of  G-  In  order  to  find  the  Fourier  transform  of  G 
we  need  to  introduce  the  hydrodynamical  equations  for  a  pressure  point 
moving  at  constant  velocity  across  the  water  surface.  The  hydrody¬ 
namical  problem  is  classical  (e.g. ,  Lamb,  1879)  and  has  been  treated 
by  a  number  of  authors.  Most  of  the  classical  treatments  are  con¬ 
cerned  with  the  form  of  Q.  However,  here  we  are  only  interested  in 
its  transform  Z  which  is  somewhat  easier  to  derive. 
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We  assume  the  pressure  point  moves  with  velocity  v  m  the 
negative  x  direction  and  at  the  instant  of  observation  it  is  located  at 
the  origin.  Assuming  irrotational  flow  we  may  introduce  a  velocity 

A 

potential  CP  and  the  equations  of  hydrodynamics  become:  the  equation 
of  incompressibility 

V2  $  =  0  ,  z  <  0  (21) 

the  equation  of  motion 

A 

V  !~  =  g  C  +  ~P(x,y)  ,  Z  =  0  (22) 

and  the  condition  that  the  normal  component  of  velocity  be  zero  at  a 
fixed  boundary 

+  v  =  0  ,  z  =  0  (23) 

o  z  ox 

In  these  equations,  p  represents  the  density  of  the  water,  P  =  PQ6(x,y) 
is  the  magnitude  of  the  pressure  and  g  the  acceleration  of  gravity.  It 
should  be  pointed  out  that  (22)  and  (23)  are  linearized  approximations 
and  require  cp  and  C  be  small.  The  origin  of  the  xyz  axes  is  moving 
in  the  negative  x  direction  with  a  uniform  velocity  v.  Now  we  take  the 
two-dimensional  Fourier  transform  over  x  and  y;  i.e., 

$  (S.Tl.z)  =  J  dx  j  dy  £{x,y,z)  e  +  ^  (24) 

-  CD  -  CO 

together  with  the  Fourier  transform  of  G  given  in  (19)  and  equations 
(21),  (22)  and  (23)  become 


-(?2  +  Tl3)  *•  +  -r-jf  =  0  , 

0  z 

z  <  0 

(25) 

A  p_ 

i§v§  =  gZ  +  , 

z  =  0 

(26) 

A 

-N  X 

0  i  -si 

- —  =  -  iv?  Z 

0  z 

N 

II 

O 

(27) 
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Eliminating  Z  from  (26)  and  (27)  gives 


ifvf  i  =  ±* . 

g  S  z  g  p  »  z 

The  solution  of  (25)  is 

i  _  $0  e/ ?2  +  ri2  z 

Substituting  (29)  into  (28)  and  solving  for  ij  gives 


-i  vP 


$0  = 


gp 


_  -J  l2  +  Tl2 


F rom  (27)  we  can  solve  for  Z 


Z(S,ri)  = 


(l2  +  r\3)z 


gP  "(?2  +T13)*  ' 


(28) 


(29) 


(30) 


(31) 


The  water  wave  number  is  kw  =  g/v2  and  substituting  k  into  (31) 
yields  the  desired  result 


Z(§.T1) 


[>-kw  <5= 


(32) 


the  Fourier  transform  of  £(x, y).  The  quantity  (PW  gp)  has  dimen¬ 
sions  (length)2  and  Z  has  dimensions  (length)3  .  Substituting  (32) 
into  (18)  and  using  the  acoustical  interpretation  for  §  and  T|  given 
in  (20)  gives 
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(33) 


the  scattered  acoustic  field.  The  scattered  fiold  has  dimensions 
(length)  1  which  is  the  same  as  the  incident  field. 

We  now  turn  to  an  interpretation  of  (33).  The  denominator 
vanishes  when 

2  k  cos2  3  =  (34) 

The  resulting  singularity  perhaps  requires  further  analysis,  but  it  is 
clear  that  it  represents  a  "resonance".  For  the  moment  we  will  think 
of  and  Ti  as  components  of  the  water  wave  number  and  not  the 
components  of  the  acoustic  wave  number  as  defined  in  (20).  This 
interpretation  is  motivated  by  noting  that  the  denominator  in  (32)  is 
identical  to  that  in  (33)  and  the  result  in  (32)  has  nothing  to  do  with  the 
acoustic  problem.  It  is  the  Fourier  transform  of  the  water  surface. 
The  water  wave  phase  is 

t  =  x§  +  y  T)  (35) 

and  the  stationary  water  phase  points  are  solutions  to 


dl|> 


-  (-S-) 


(36) 


At  resonance,  from  (32),  §  and  T)  are  related  as 


§2  =  k  (52  +  T]2)2 
w 


(37) 


Solving  (35),  (36)  and  (37)  for  x  and  y  in  terms  of  a  parameter 

u  =  §  /  k  gives 
w 
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! 


x  =  a  L^3  -!) 
u3 


y  =  +  a 


(u3  -1) 
if3 


(38) 

(39) 


with  a  =  ^  /  kw.  These  two  equations  represent  the  surfaces  of  constant 
water  wave  phase.  One  water  wave  crest  is  shown  in  figure  2  where 
a  =  Xw.  To  generate  the  next  crest  a  =  2^w  and  in  general  the  pth 
crest  is  generated  by  a  =  p\  w>  The  cuspid'll  lines  forming  the  boundary 
of  the  sector  occupied  by  the  wake  may  be  determined  from  (38)  and 
(39)  by  defining 


tana 


=  Z  =  (u3  -*) 


x 


2u3  -  1 


(40) 


and  finding  the  point  where 


5  (tana) 
du  =  0 


(41) 


Substituting  (40)  into  (41)  gives  the  point 

u  =  i  J1  (42; 

Substituting  (42)  back  into  (40)  gives  the  angular  sector  defining  the 
wake  (Kelvin,  1887) 


a  =  +  tan 
c  — 


1  fey?)*1’9’28' 


(43) 


x  and 


This  cuspidal  angle  is  shown  in  figure  2.  For  1  £  u  £  J  — 

I  y|  lie  on  the  portion  of  the  curve  starting  from  the  point  x=  1,  y  =  0 
and  going  to  the  cusp  (transverse  component).  For  u  £  ,  x  and 

I  y|  he  on  the  portion  of  the  curve  between  the  cusp  and  the  origin 
(lengthwise  component). 
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The  tangent  to  the  curve  in  figure  2  is  given  by 


dy  _  dy/  du  _  I 
dx  dx/  du  +  (u3  -1)2 


(44) 


The  tangent  to  the  cusps  in  figure  2  are  given  by  (44)  where  u2  =  3/2 


aT  =  tan"1  J  2  ^  +  54*  44# 


(45) 


This  angle  is  shown  in  figure  2.  The  angles  and  aT  satisfy 

a  =  2a  -  n/  2. 

c  T 


Finally  ,  we  return  to  the  interpretation  of  (3  3)  from  the  scat¬ 
tering  viewpoint.  Consider  our  acoustic  wave  incident  at  an  angle  0 
with  respect  to  the  negative  x-axis  as  shown  in  figure  3.  A  tangent 
vector  to  the  water  wave  phase  surfaces  in  figure  3  (shown  as  a  unit 
vector  e^,  ,  in  figure  3)  is 


u^  1 


(46) 


The  incident  acoustic  wave  will  be  pci pendicular  to  the  water  wave 
phase  surfaces  giving  maximum  reflection  in  the  backscatter  direction 
(resonance)  when 


T  ' 


=  0 


with 


e.  =  e  cos  9  -  e  sin  9 
—  k  — x  — y 


(47) 


(48) 


Substituting  (46)  and  (48)  into  (47)  gives 

I 

tan  9  =  (u2  -1)  2  (49) 

It  is  not  surprising  that  the  result  in  (49)  is  the  reciprocal  of  that  in 
(44).  The  surfaces  of  constant  wa.ter  wave  phase  can  be  thought  of  as 
representing  a  diffraction  grating.  When  twice  the  path  length  d  (to 
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go  in  and  out)  is  equal  to  the  incident  acoustic  wave  length  A.  ,  neighbour¬ 
ing  water  wave  lines  in  the  grating  will  interfere  constructively  pro¬ 
ducing  maximum  backscatter  (resonance).  Twice  the  distance  d  in 
figure  3  is 

2d  =  2  R  •  ek  (50) 

where 


R  =  xe  +  y  e 

-  -x  -y 


2tt 


( 2u2 - 1 ) 
uJ 


+  e 

-y 


Substituting  (48)  and  (51)  into  (50)  yields 


„  ,  V  w  ,  n  (2u2  -1)  .  a  (u2-l)2 

2d  =  ~  lcos0  — —  - sm0  — S 


This  distance  remains  the  same  for  any  two  adjacent  water  wave 
surfaces  in  the  grating.  Substituting  for  u  in  (52)  in  terms  of  0  from 
(49)  gives  ^ 

2d  =  cos20  (53) 

The  phase  difference  'll  between  two  adjacent  water  wave  crests  is 
2tt  and  (53)  becomes 

2d  =  Z\  cos20  (54) 

w 

or 

A.  =  2  K w  cos20  (55) 

which  is  identical  with  (34). 
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3:  DISCUSSION 


An  analysis  for  scattering  of  scalar  waves  in  air  by  water 
waves  is  developed.  In  the  analysis  we  have  made  the  assumption  that 
the  wake  extends  to  infinity.  In  cases  where  noise  or  clutter  is  pre¬ 
sent  the  wake  would  attenuate  and  would  not  continue  indefinitely  behind 
the  pressure  point.  The  effect  of  clutter  or  dissipation  may  be  intro¬ 
duced  into  this  analysis;  however,  the  condition  for  resonance  given 
in  (55)  remains  unaltered.  The  dissipation  factor  that  is  introduced 
into  the  hydrodynamical  equations  (21) -(23)  eliminates  the  standing 
wave  that  extend  an  infinite  distance  in  back  of  the  pressure  point.  By 
considering  the  scattered  field  as  a  generalized  function  (Hufford, 

1972)  and  taking  the  limit  as  the  dissipation  tends  to  zero  we  obtain  the 
result  given  in  (55). 

The  result  in  (55)  is  interesting  from  another  viewpoint.  Kurss 
and  Crombie  (1971)  derived  the  following  resonance  condition 

\  =  2  cos(9-0tc)  (56) 

Their  analysis  was  based  on  examining  the  behavior  of  the  scattered 
field  only  in  the  vicinity  of  the  cuspidal  lines.  From  figure  2  we  see 
that  the  cuspidal  region  is  only  a  small  portion  of  the  wake.  An  excel¬ 
lent  photograph  showing  the  intensity  of  the  wake  near  the  track  of  the 
pressure  point  is  given  in  Feynman  et  al.  (1963).  This  photograph 
shows  the  intensity  near  the  track  to  be  nearly  as  great  as  the  intensity 
near  the  cuspidal  lines.  Numerically  generated  wakes  (Tuck  et  al.  , 
1971)  also  show  a  strong  component  of  wake  along  the  track.  Thus,  it 
is  not  surprising  that  the  resonance  relation  based  on  the  entire  wake 
is  different  from  the  resonance  condition  based  upon  the  cuspidal 
component  of  the  wake. 


-13- 


4.  REFERENCES 

Feynman,  R.  P.  ,  Leighton,  R.  B.  and  Sands,  M.  (1963),  The  Feynman 
Lectures  on  Physics,  Vol.  1,  Addison-Wesley  Publishing  Co. 
(New  York),  p.  51-8. 

Havelock,  T.  H.  (1918),  Wave  resistance:  some  cases  of  three- 
dimensional  fluid  motion,  Proc.  of  the  Royal  Society,  A.  , 

Vol.  95,  Nov.  27,  p.  147,  equations  (1)  and  (2). 

Hogner,  E.  (1923),  A  contribution  to  the  theory  of  ship  waves,  Ark. 
Mat.  Astr.  Fys.,  17,  No.  12,  pp.  1-50. 

Hufford,  G.  A.  (1972),  Scattering  by  a  perturbed  plane  and  generalized 
functions,  unpublished  notes. 

Lord  Kelvin  (1887),  On  ship  waves,  Proc.  of  the  Royal  Scoeity, 

Vol.  42,  p.  80. 

Kurss,  H.,  and  Crombie,  D.  D.  (1971),  private  communication. 

Sir  Horace  Lamb  (1879)*  Hydrodynamics,  Dover  Publications, 
pp.  384-440. 

Peters,  A.  S.  (1949),  A  new  treatment  of  the  ship  wave  problem. 
Comm.  Pure  Appl.  Math.,  Vol.  2,  pp.  123-148. 

Sommerfeld,  A.  (1950),  Mechanics  of  deformable  bodies,  Academic 
Press  (New  York),  Ch.  5. 

Tuck,  E,  O.,  Collins,  J.  I.  and  Wells,  W.  H.  (1971),  On  ship  wave 
patterns  and  their  spectra.  Jour,  of  Ship  Research,  Vol.  15, 

No.  1,  pp.  11-21. 

Ursell,  F.  (1960),  On  Kelvin' s  ship-wave  pattern,  Jour,  of  Fluid 
Mech.  ,  Vol.  8,  pp.  418-431. 


-14- 


5.  ACKNOWLEDGEMENT 


The  authors  wish  to  thank  Dr.  Donald  E.  Barrick,  Dr.  James 
R.  Wait  and  Mr.  D.  D.  Crombie  for  their  many  constructive  sugges¬ 
tions  during  this  investigation.  The  authors  also  thank  Mr.  Alvin  Van 
Every  for  suggesting  the  problem. 


P,  Source  and 
Observer 


z  =  C  U,y) 

0,  Integration  Point 

► 


Figure  1.  Geometry  for  scattering  problem. 
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